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In  this  article,  we  present  a  quantitative  stiffness  imaging  technique  and  demonstrate  its  use  to 
directly  map  the  dynamic  mechanical  properties  of  materials  with  nanometer-scale  lateral 
resolution.  For  the  experiments,  we  use  a  “hybrid”  nanoindenter,  coupling  depth-sensing 
nanoindentation  with  scanning  probe  imaging  capabilities.  Force  modulation  electronics  have  been 
added,  enhancing  instrument  sensitivity  and  enabling  measurements  of  time  dependent  materials 
properties  (e.g.,  loss  modulus  and  damping  coefficient)  not  readily  obtained  with  quasi-static 
indentation  techniques.  Tip-sample  interaction  stiffness  images  are  acquired  by  superimposing  a 
sinusoidal  force  (~1  fiN)  onto  the  quasi-static  imaging  force  (1.5-2  /xN),  and  recording  the 
displacement  amplitude  and  phase  as  the  surface  is  scanned.  Combining  a  dynamic  model  of  the 
indenter  (having  known  mass,  damping  coefficient,  spring  stiffness,  resonance  frequency,  and 
modulation  frequency)  with  the  response  of  the  tip-surface  interaction,  creates  maps  of  complex 
stiffness.  We  demonstrate  the  use  of  this  approach  to  obtain  quantitative  storage  and  loss  stiffness 
images  of  a  fiber-epoxy  composite,  as  well  as  directly  determine  the  loss  and  storage  moduli  from 
the  images  using  Hertzian  contact  mechanics.  Moduli  differences  as  small  as  20%  were  resolved  in 
the  images  at  loads  two  orders  of  magnitude  lower  than  with  indentation,  and  were  consistent  with 
measurements  made  using  conventional  quasi-static  depth-sensing  indentation  techniques.  ©  2001 
American  Institute  of  Physics.  [DOI:  10.1063/1.1380218] 


I.  INTRODUCTION 

The  objective  of  the  present  work  is  to  map  mechanical 
properties  quantitatively  over  an  area  on  a  sample  surface.  In 
the  mechanics  community,  a  traditional  approach  involves 
the  use  of  indentation  arrays  to  construct  hardness  and  modu¬ 
lus  maps.1'2  However,  for  quantitative  and  meaningful  mea¬ 
surements,  the  indents  have  to  be  finitely  spaced  and  thou¬ 
sands  of  indents  have  to  be  performed  to  obtain  even  a 
rudimentary  map;  this  is  time  consuming  and  lateral  resolu¬ 
tion  is  limited  by  the  indent  spacing. 

Mechanical  property  mapping  can  be  significantly  en¬ 
hanced  by  the  addition  of  tip-sample  rastering  technology 
(through  piezoelectric  scanners),  enabling  much  finer  tip- 
sample  positioning  below  the  optical  limit.  With  cantilever- 
based  atomic  force  microscopy  (AFM),  maps  of  adhesion3  5 
and  elastic  properties6,7  can  be  obtained  through  multiple  in¬ 
dividual  force-distance  curves.  Resolution  is  again  limited 
by  contact  area  if  the  tip-sample  interaction  is  elastic,  or  by 
appropriate  spacing  if  the  deformation  is  plastic.  Quantitative 
analyses  of  AFM  data  are  complicated  by  problems  relating 
to  tip  deformation  and  geometry,  piezo  hysteresis  and  creep, 
as  well  as  an  off-normal  tip  approach.*-10 

More  recently,  numerous  AFM  imaging  techniques  em¬ 
ploying  modulation  have  been  developed.  In  these  tech- 
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niques,  the  changes  in  amplitude  or  phase  response  of  an 
oscillating  tip  or  sample  as  the  tip  is  scanned  over  a  hetero¬ 
geneous  sample  can  be  used  to  map  materials  property  varia¬ 
tions.  Most  of  these  techniques  utilize  displacement  modula¬ 
tion,  where  either  the  tip  holder  or  sample  is  oscillated,  and 
the  tip-sample  response  (amplitude  or  phase  shift)  used  to 
provide  an  image.11-17  In  these  techniques,  the  tip  can  be 
either  in  continuous  or  intermittent  contact  with  the  sample. 

While  these  configurations  are  powerful  imaging  tools, 
quantitative  interpretation  of  the  data  requires  knowledge  of 
often  complicated  system  dynamics.18  Coupled  with  the 
complex  nature  of  the  tip-sample  interaction  (especially 
with  intermittent  contact  modes),  the  resulting  data  are  diffi¬ 
cult  to  model  and  often  misinterpreted.  For  example,  phase 
shifts  observed  in  one  form  of  intermittent  contact  mode  im¬ 
aging  (“tapping  mode”)  are  necessarily  due  to  energy  dissi¬ 
pation  (damping  and  adhesion);19,20  These  dissipative  pro¬ 
cesses  are  often  ignored  and  the  phase  response  incorrectly 
ascribed  purely  to  elasticity  variations.  Hence,  while  imaging 
“elasticity”  is  possible,  quantification  is  nontrivial.  Further¬ 
more,  the  inherent  problems  associated  with  vibrating  canti¬ 
levers  (e.g.,  bistabilities21  and  multiple  resonances),  un¬ 
known  electronics  corrections,  and  feedback  issues  make 
quantitative  interpretation  of  image  data  extremely 
difficult.22,23 
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Like  normal  modulation  techniques  (e.g.,  perpendicular 
to  sample  normal),  lateral  modulation  of  the  sample  in  the 
AFM  can  be  used  to  map  sample  inhomogeneities  like 
defects12  or  friction  contrast,24  as  well  as  glass  transition 
temperatures.25  If  great  care  is  taken,  quantitative  stiffness 
images  can  be  obtained;26  however,  potential  frictional  con¬ 
tributions  due  to  slip,  as  well  as  cantilever  issues  (e.g., 
matching  cantilever  and  sample  stiffnesses,  and  lateral  force 
calibration  difficulties)  limit  the  technique. 

Force  modulation  methods,  where  a  force  (rather  than 
displacement)  is  applied  directly  to  the  tip,  have  been  used  in 
indentation  to  enhance  sensitivity.  If  the  dynamic  response  of 
the  indenter  is  well  calibrated  and  modeled,  dynamic  me¬ 
chanical  properties  (e.g.,  storage  and  loss  moduli)  can  be 
measured.  ’  Force  modulation  is  also  possible  in  AFM,  ’ 
which  gives  far  greater  sample  positioning  and  added  imag¬ 
ing  capabilities  compared  to  conventional  nanoindenters. 
However,  the  small  forces  applied  in  AFM  limit  the  range  of 
samples  that  can  be  studied. 

The  sample  positioning  limitations  of  conventional  in¬ 
dentation  techniques  can  be  overcome  through  hybrid  indent- 
ers,  coupling  depth-sensing  indentation  with  piezo  scanning 
capabilities.31'32  With  a  hybrid  indenter,  AFM-like  intermit¬ 
tent  contact  imaging  using  force  modulation  has  been 
demonstrated.33  More  recently,  we  have  shown  that  with 
contact  force  modulation  using  a  hybrid  nanoindenter,  it  is 
possible  to  quantitatively  map  the  surface  stiffness.34'35 
These  stiffness  images  are  produced  by  monitoring  the  tip 
displacement  amplitude  and  phase  responses  to  an  applied 
sinusoidal  force,  and  are  properly  interpreted  as  a  complex 
quantity  containing  both  elastic  (storage)  and  dissipative 
(loss)  components. 

In  this  article,  we  present  an  imaging  technique  that  we 
use  to  map  the  tip-sample  mechanical  response  quantita¬ 
tively  using  the  hybrid  indenter.  The  advantage  of  this  ap¬ 
proach  is  that  not  only  can  it  be  used  to  map  the  stiffness 
contrast  but  also  to  map  the  elastic  modulus  quantitatively, 
for  both  storage  and  loss  components.  The  storage  and  loss 
moduli  determined  from  the  images  are  compared  to  the 
measurements  made  with  conventional  depth-sensing 
nanoindentation  techniques.  Issues  that  must  be  addressed  in 
order  to  map  the  mechanical  properties  quantitatively,  such 
as  instrument  dynamics,  feedback  modulation  of  the  piezo, 
indenter  tip  shape  calibration,  as  well  as  contributions  of 
adhesion  and  sliding,  are  discussed. 

II.  EXPERIMENT 


used  to  measure  the  amplitude  and  phase  of  the  displacement 
signal.  Finally,  to  process  the  data,  appropriate  contact  me¬ 
chanical  models  and  knowledge  of  tip  geometry  are  neces¬ 
sary.  Our  approach  to  solving  these  problems  is  described 
next.  We  demonstrate  measurement  of  storage  and  loss  con¬ 
tact  stiffness  during  scanning  in  continuous  contact  with  a 
sample,  resulting  in  images  that  map  storage  and  loss  com¬ 
ponents  of  the  elastic  modulus. 

The  instrument  configuration  used  in  these  experiments 
consists  of  a  hybrid  indenter  (Triboscope-Hysitron,  Inc., 
Minneapolis,  MN)  on  an  AFM  base  (Multimode  with  JV 
scanner  and  signal  access  module  -  Digital  Instruments, 
Santa  Barbara,  CA),  with  added  ac  force  modulation  elec¬ 
tronics  described  previously.36  Briefly,  the  indenter  trans¬ 
ducer  consists  of  a  capacitive  core  having  two  fixed  outer 
plates  and  a  spring  suspended  center  plate  to  which  the  in¬ 
denter  tip  is  attached.  The  center  plate  is  actuated  by  attract¬ 
ing  the  center  plate  towards  the  fixed  lower  plate  by  electro¬ 
static  means,  and  the  resulting  displacement  detected  by 
capacitive  sensing.  Coupling  of  the  transducer  to  the  AFM 
scanner  allows  both  depth-sensing  indentation  (with  the  pi¬ 
ezo  feedback  turned  off)  and  topographic  imaging  (using  the 
piezo  feedback).  Our  method  of  dynamic  mechanical  prop¬ 
erty  mapping  is  described  next. 

A.  Dynamic  indentation:  Mechanical  and  electronic 
contributions 

To  deconvolute  the  actual  dynamic  response  of  the  con¬ 
tact  in  a  dynamic  measurement,  the  contributions  from  the 
various  electronics  (e.g.,  filters,  etc.)  and  mechanical  ele¬ 
ments  must  be  known  and  modeled  since  they  will  be  depen¬ 
dent  on  the  frequency  of  the  measurement.  This  problem  can 
be  solved  if  the  transfer  functions  for  the  electronics  and 
mechanical  elements  of  the  transducer  are  known  or  cali¬ 
brated. 

In  our  previous  work,  we  have  shown  that  the  dynamic 
mechanical  response  of  the  transducer  when  it  is  in  contact 
with  the  sample  can  be  modeled  using  two  Kelvin-Voigt 
mechanical  equivalents.36  For  contact  stiffness,  much  less 
than  the  machine  stiffness  (at  least  two  orders  of  magnitude), 
the  machine  can  be  assumed  to  be  rigid  and  the  model  solved 
analytically.27  In  this  case,  the  amplitude  of  the  displacement 
signal  is 


\J(k  —  m  a>2)2  +  \(C  j+  C  s)  a>]2 
and  the  phase  shift  between  force  and  displacement  is 


One  of  the  barriers  preventing  quantitative  mechanical 
property  imaging  in  AFM  or  other  scanning  techniques  is 
that  the  measured  response  is  the  convolution  of  the  tip- 
sample  interaction  and  the  measuring  instrument.  A  number 
of  factors  must  be  accounted  for  to  interpret  this  response. 
For  dynamic  measurements,  the  system  mechanics  (mass, 
spring  constant,  and  damping  loss)  and  instrument  compli¬ 
ance  must  be  modeled  and  measured,  and  system  electronics 
corrections  must  be  made.  Acquiring  this  information  while 
scanning  adds  the  need  to  control  the  system  feedback  re¬ 
sponse,  scan  rates,  and  time  constant  of  the  lock-in  amplifier 


(/>  =  tan 


(C,  +  Cs)w 
k  —  mur 


(2) 


where  X{)  is  the  amplitude  of  the  displacement  oscillation 
and  cf>  is  the  phase  shift  of  the  displacement  with  respect  to 
the  exciting  force,  F0  is  the  amplitude  of  the  sinusoidal 
force,  m  is  the  indenter  mass,  oj  is  the  frequency  in  radians, 
C,  is  the  damping  coefficient  of  the  air  gap  in  the  capacitive 
displacement  sensor,  and  Cs  is  the  damping  coefficient  of  the 
contact.  The  combined  stiffness,  k,  is  given  by 

k  =  Ks  +  Kj ,  (3) 
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where  Ks  is  the  tip-sample  contact  stiffness  and  K ,  is  the 
spring  constant  of  the  leaf  springs  that  hold  the  indenter 
shaft. 

The  indenter  spring  constant  Kj ,  mass  m  and  damping 
coefficient  of  the  air  gap  C,  can  be  obtained  by  fitting  the 
amplitude  and  phase  response  of  the  instrument  (between  10 
and  300  Hz)  to  Eqs.  (1)  and  (2)  when  it  is  not  in  contact  with 
a  sample.  However,  for  an  accurate  calibration,  the  contribu¬ 
tion  due  to  the  electronics  must  be  determined.  In  our  instru¬ 
ment,  the  electronic  contributions  to  the  dynamic  response 
are  due  to  low-pass  filters  in  the  force  actuation  and  displace¬ 
ment  sensing  electronics.  If  the  transfer  functions  of  the  elec¬ 
tronics  are  known,  the  electronic  contribution  for  the  dy¬ 
namic  response  can  be  corrected  (see  Appendix).  Figure  1(a) 
shows  the  combined  dynamic  compliance  of  the  instrument 
(i.e.,  ratio  of  the  displacement  amplitude  X0  to  the  force 
amplitude  F0  with  electronic  contribution  directly  measured 
from  the  lock-in  amplifier)  when  it  is  free  of  contact  with  the 
sample  without  any  correction  for  electronics.  If  we  attempt 
to  fit  the  response  to  Eqs.  (1)  and  (2),  it  results  in  a  poor  fit 
compared  to  the  mechanical  constants  obtained  from  quasi¬ 
static  (dc)  measurements. 

Figure  1  (b)  shows  the  dynamic  compliance  after  correct¬ 
ing  for  the  electronic  contribution  using  Eqs.  (Al),  (A3),  and 
(A5)  from  the  appendix.  The  fit  for  the  indenter  constants 
was  Kj  (133  N/m),  m  (276  mg),  and  C(-  (0.0208  Ns/m);  the 
values  for  Kj  and  m  are  in  good  agreement  with  quasi-static 
measurements.  Once  the  contributions  from  the  electronics 
are  corrected,  the  measured  system  response  correlated  well 
with  the  Kelvin- Voight  mechanical  equivalent  model.  Figure 
1(c)  shows  the  difference  between  the  measured  phase  shift 
response  cj)  with  and  without  electronic  correction.  For  the 
electronic  correction,  the  total  phase  shift  due  to  electronics 
[Eqs.  (A2),  (A4),  and  (A6)]  is  subtracted  from  the  phase  shift 
measured  by  the  lock-in  amplifier. 

The  machine  compliance  was  measured  (using  the 
method  of  Oliver  and  Pharr)28  and  found  to  be  25  m//xN.  If 
all  the  aforementioned  constants  are  known  or  accurately 
calibrated,  the  contact  stiffness  Ks  and  contact  damping  co¬ 
efficient  Cs  can  be  calculated  using  Eqs.  (l)-(3).  In  an  in¬ 
dentation  experiment,  the  contact  stiffness  is  proportional  to 
the  projected  contact  area.  A, 

K=2E*J~,  (4) 

"  7 T 

where  E*  is  the  reduced  modulus  (combining  tip  and 
sample).  For  materials  which  can  damp  the  superimposed 
dynamic  force,  K's  is  the  storage  component  of  the  contact 
stiffness  in  phase  with  the  applied  force,  and 

K"s  =  u>Cs,  (5) 

where  K"  is  the  loss  component  of  the  contact  stiffness  out 
of  phase  with  the  applied  force.  Storage  and  loss  moduli  for 
polymer  materials  obtained  by  this  technique37  39  compare 
favorably  with  bulk  values  determined  through  conventional 
dynamic  testing. 

The  lock-in  amplifier  also  contributes  a  time  constant  to 
the  measurement.  The  settings  are  usually  a  compromise  be¬ 
tween  tolerable  noise  and  speed  of  data  acquisition.  For  ex- 


Frequency,  radians 


FIG.  1.  Plots  of  dynamic  compliance  (amplitude  ratio  of  X0/F0)  of  the 
indenter  system  including  (a)  electronic  contribution  and  (b)  after  correcting 
for  the  electronic  contribution  are  shown.  The  solid  line  is  the  fit  for  the 
Kelvin-Voight  model.  The  corrected  measurement  follows  the  model  pre¬ 
diction.  (c)  The  phase  shift  as  a  function  of  frequency  with  and  without 
electronic  contribution. 


ample,  in  a  static  measurement  (i.e.,  at  fixed  load),  one  could 
optimize  for  the  low  noise  in  the  measurement.  However, 
both  indentation  and  scanning  operations  require  knowledge 
of  the  dynamic  response  of  the  indenter  to  the  small  sinu¬ 
soidal  force  while  moving  the  sample  vertically  or  laterally 
with  respect  to  the  indenter  tip.  The  frequency  of  the  mea¬ 
surement  and  rate  of  sample  movement  will  determine  ac¬ 
ceptable  time  constant  settings.  If  the  time  constant  setting 
on  the  lock  in  is  too  high,  the  data  will  be  averaged  over 
longer  time  scales  than  the  apparent  resolution  of  the  mea¬ 
surement,  distorting  the  response.  In  this  case,  a  purely  elas- 
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FIG.  2.  The  amplitude  of  the  indenter  tip  response  as  a  function  of  fre¬ 
quency  for  three  different  conditions:  •  when  the  tip  is  not  in  contact  with 
the  surface,  ■  when  the  tip  is  in  contact  with  the  surface  (~1  /xN  load) 
under  feedback  control,  and  ♦  when  the  tip  is  in  contact  with  the  sample 
(~1  /uN  load)  without  feedback  are  shown. 


tic  deformation  may  appear  plastic,  while  an  image  could  be 
blurred  spatially.  Care  must  be  taken  to  ensure  that  the  speed 
of  the  experiment  is  not  greater  than  the  rate  at  which  the 
lock  in  can  detect  variation  in  the  response  of  the  instrument. 


B.  Feedback  issues 

The  use  of  feedback  electronics  during  imaging  presents 
a  significant  issue  influencing  quantitative  mechanical  mea¬ 
surements.  The  resonance  frequency  of  our  instrument  (110 
Hz)  and  the  useable  dynamic  range  for  force  modulation 
(10-250  Hz)  are  below  the  frequency  of  the  feedback  con¬ 
troller  (in  the  kHz  range)  for  constant  load  scanning.  There¬ 
fore,  when  the  tip  is  modulated,  the  feedback  control  will 
force  the  piezo  scanner  holding  the  sample  to  counter  the 
vertical  oscillation  of  the  tip  to  keep  constant  deflection,  an 
undesirable  response.  Figure  2  shows  the  amplitude  of  the 
indenter  tip  as  a  function  of  frequency  for  three  different 
conditions:  •  when  the  tip  is  not  in  contact  with  the  surface, 
■  when  the  tip  is  in  contact  with  the  surface  (~1  /xN  load) 
with  feedback  control,  and  ♦  when  the  tip  is  in  contact  with 
the  sample  (~1  /xN  load)  without  feedback.  When  the  tip  is 
in  contact,  the  resonance  frequency  increases  (above  the  free 
resonance  frequency  of  110  Hz).  At  low  frequencies,  the  tip 
displacement  amplitude  is  lower  when  the  tip  is  in  feedback 
than  without  feedback.  The  result  is  an  apparent  “stiffening” 
of  the  tip-sample  interaction,  increasing  the  measured  con¬ 
tact  stiffness;  this  would  be  expected  since  the  feedback  is 
counteracting  the  tip  oscillation.  The  higher  amplitudes  ob¬ 
served  at  a  higher  frequency  with  feedback  on  is  likely  due 
to  the  effects  of  the  overall  system  damping.  We  also  observe 
a  change  in  phase  shift  under  feedback  control. 

To  overcome  these  problems,  we  have  implemented  an 
interleaved  scanning  technique  provided  with  the  AFM  scan¬ 
ning  software.  During  the  first  pass,  the  feedback  control  is 
active  and  the  height  information  is  stored.  Then,  during  the 
next  pass,  the  feedback  is  turned  off  and  the  height  data 
collected  from  the  first  pass  is  used  to  position  the  Z  piezo  to 
follow  the  surface  topography.  At  the  same  time,  the  tip  dis¬ 
placement  amplitude  and  phase  data  from  the  lock-in  ampli¬ 
fier  and  the  contact  force  data  from  the  displacement  trans¬ 


ducer  are  collected.  In  this  way,  the  phase  and  amplitude 
response  of  the  contact  to  the  sinusoidal  force  can  be  accu¬ 
rately  measured. 

C.  Sample  preparation 

A  composite  made  of  carbon  fibers  (nominally  6  /xm  in 
diameter)  embedded  in  an  epoxy  matrix  was  used  for  the 
experiments.  A  2  mm  thick  sample  was  cut  using  a  diamond 
saw  and  mechanically  polished;  the  final  polish  was  per¬ 
formed  using  0.25  /xm  diamond  paste.  The  sample  was  then 
mounted  onto  the  AFM  piezo  scanner  inside  a  dry  box 
purged  with  dry  nitrogen  (RH<2%).  To  minimize  thermal 
drift,  the  sample  was  allowed  to  equilibrate  with  the  tem¬ 
perature  of  the  dry  box  for  at  least  1  h. 

D.  Scanning  conditions  for  stiffness  images 

Diamond  indenters  with  tip  radii  between  250-600  nm 
were  used  for  the  experiments.  The  indenter  tips  (tip  radius  R 
and  area  function)  were  calibrated  using  both  standard  inden¬ 
tation  techniques28  with  quartz  and  blind  tip 
reconstruction40'41  with  a  tip  characterizer  (TGTOl-Silicon- 
MDT,  Moscow,  Russia).  The  sample  was  then  brought  into 
contact  with  the  diamond  tip  at  a  constant  (dc)  contact  force 
of  1.5-2  /xN  using  AFM  feedback  control.  To  obtain  the 
stiffness  image,  the  indenter  tip  was  modulated  with  a  sinu¬ 
soidal  force  of  ~1  /xN  at  200  Hz,  which  resulted  in  a  tip 
displacement  amplitude  between  0.5-2  nm.  Interleaved 
scanning  was  performed  at  a  scan  rate  of  0.1  Hz  over  100 
/xm2  and  512X512  pixel  resolution  for  all  the  experiments. 
The  lock-in  time  constant  was  set  to  3  ms.  Data  were  ac¬ 
quired  and  processed  to  obtain  contact  force,  sinusoidal  dis¬ 
placement  amplitude  and  phase  shift.  Electronics  correction 
factors  (at  200  Hz)  were  0.594  for  force  [see  Eq.  (Al)], 
1.163  for  displacement  (1/[(A3)*(A5)]  and  66°  for  the  phase 
shift.  The  data  were  further  processed  using  image  manipu¬ 
lation  software  (Transform,  Research  Systems,  Boulder,  CO) 
to  obtain  the  damping  coefficient  as  well  as  storage  and  loss 
stiffness  images  using  Eqs.  (l)-(3). 

E.  Indentation  experiments 

Both  standard  force-displacement28  (to  a  maximum  load 
of  200  /xN)  and  force-modulation  methods27  were  used  to 
evaluate  moduli  of  the  various  regions  of  the  sample.  Stan¬ 
dard  calibration  techniques  were  employed  to  model  the  area 
function  of  the  Berkovich  tip  for  indents  to  depths  greater 
than  the  tip  defect.28  The  force-modulation  amplitude  for  in¬ 
dents  ranged  between  1-3  /xN  at  105  Hz,  with  a  lock-in  time 
constant  of  100  ms  and  indentation  rate  of  2-20  /xN/s. 
Force-distance  curves42  were  used  to  evaluate  the  adhesion 
of  the  indenter  tip  to  the  sample. 

III.  RESULTS  AND  DISCUSSION 
A.  Stiffness  imaging  and  interpretation 

Figure  3(a)  shows  the  contact  force  image  obtained  by 
monitoring  the  indenter  tip  displacement  using  the  displace¬ 
ment  transducer  and  the  known  spring  constant  of  the  trans¬ 
ducer.  Although  the  contact  force  was  kept  constant  around 
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FIG.  3.  Images  of  a  carbon  fiber-epoxy  composite,  scan  size  is  10  /imXlO  /Am,  for  (a)  the  contact  force  (in  /a N),  (b)  the  ratio  of  the  force  modulation 
amplitude  to  the  displacement  amplitude  (in  N/m  and  corrected  for  electronics,  see  the  Appendix)  corresponding  to  the  contact  force  image,  and  (c)  the  true 
phase  shift  (in  degrees  and  corrected  for  electronics)  between  the  oscillating  force  and  the  displacement  during  scanning  are  shown. 


1.6  yuN,  a  small  variation  (100  nN)  can  be  seen  as  a  result  of 
feedback  error,  giving  slight  contrast  to  the  image.  Figure 
3(b)  shows  the  ratio  of  the  sinusoidal  force  amplitude  F0  to 
the  displacement  amplitude  X0  (corrected  for  electronics,  see 
the  Appendix)  corresponding  to  the  constant  contact  force 
imaging.  The  measured  amplitude  ratio  is  due  primarily  to 
the  mechanical  response  of  the  sample  to  the  applied  load. 
Figure  3(c)  shows  the  true  phase  shift  (corrected  for  electron¬ 
ics)  between  the  sinusoidal  force  and  the  displacement  dur¬ 
ing  scanning.  From  Eqs.  (1)  and  (2),  it  can  be  observed  that 
amplitude  ratio  and  phase  shifts  can  occur  as  a  result  of 
changes  in  the  local  damping  coefficient  and/or  contact  stiff¬ 
ness;  contrast  in  the  images  reflects  variation  of  the  com¬ 
bined  tip-sample  elastic  and  damping  response  across  the 
sample. 

The  stiffness  and  damping  response  can  be  deconvolved 
from  Figs.  3(b)  and  3(c)  using  Eqs.  (1)  and  (2).  In  this  way, 
the  images  can  be  processed  to  obtain  the  storage  and  loss 
components  of  the  contact  stiffness  using  Eqs.  (3)  and  (5). 
Figure  4  shows  the  storage  and  loss  component  of  the  con¬ 
tact  stiffness.  The  contact  stiffness  is  much  higher  for  the 
carbon  fiber  than  the  epoxy  matrix.  The  loss  components  for 
both  are  roughly  1/5  to  1/3  the  storage  components.  The 
contrast  in  both  images  is  now  shown  to  be  due  to  a  variation 
in  tip-sample  interaction  stiffness  across  the  composite 
sample,  and  is  similar  for  both  the  storage  and  loss  compo¬ 


nents.  Additionally,  a  subtle  gradient  in  contrast  at  the  center 
of  the  fiber  indicates  that  the  center  is  more  compliant  than 
the  periphery.  The  contrast  was  observed  in  the  center  of  all 
fibers,  and  could  be  due  to  porosity  or  other  nonuniformity 
often  observed  in  such  fibers.43 


B.  Materials  property  evaluation  from  images 

The  storage  and  loss  moduli  were  obtained  from  the 
stiffness  image  by  applying  an  appropriate  contact  mechan¬ 
ics  model.  For  an  elastic  contact,  in  the  absence  of  adhesion 
or  if  the  applied  force  is  much  greater  than  the  adhesion, 
Hertzian  mechanics  can  be  used  to  determine  the  contact 
radius.44  For  the  present  experiments,  force-distance  curves 
were  used  to  evaluate  the  tip-sample  adhesion.  Figure  5 
shows  a  force-distance  curve  (force  and  interaction  stiffness 
versus  sample  displacement)  for  the  carbon  fiber  region  and 
a  diamond  indenter  tip  with  a  radius  of  250  nm.  The  experi¬ 
ment  was  conducted  at  a  relative  humidity  of  —  1%.  During 
the  approach,  the  indenter  experienced  an  attractive  force  of 
—  150  nN  before  making  contact  with  the  sample.  During 
pull-off  measurements,  the  adhesion  was  150  ±30  nN.  The 
interpretation  of  the  force  and  interaction  stiffness  curves  is 
given  elsewhere.35'42  For  the  images  shown,  adhesion  from 
different  regions  of  the  sample  ranged  between  100-200  nN. 
Because  the  applied  contact  force  was  —1600  nN,  it  is  rea- 
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FIG.  4.  Images  of  a  carbon  fiber-epoxy  composite,  scan  size  is  10  /rniXlO 
/mi,  for  (a)  the  storage  and  (b)  loss  components  of  the  contact  stiffness  are 
shown. 


FIG.  5.  Force-distance  curves  showing  the  interaction  stiffness  between  a 
diamond  indenter  tip  (with  a  radius  of  250  nm)  and  carbon  fiber  are  shown. 
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FIG.  6.  (a)  Map  of  storage  modulus  (in  GPa)  and  (b)  line  scan  of  the  carbon 
fiber-epoxy  composite  are  shown.  Scan  size  is  10  //.mXlO  //m.  The  image 
contrast  is  due  to  change  in  elastic  modulus.  Line  scan  through  the  center  of 
the  image  shows  storage  modulus  gradient  at  the  center  of  the  fiber. 
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FIG.  7.  (a)  Map  of  loss  modulus  (in  GPa)  and  (b)  line  scan  of  the  carbon 
fiber-epoxy  composite  are  shown.  Scan  size  is  10  //mXlO  /mi.  Contrast  in 
the  image  is  due  to  change  in  loss  modulus.  Line  scan  through  the  center  of 
the  image  shows  loss  modulus  gradient  at  the  center  of  the  fiber. 


sonable  to  apply  Hertzian  contact  mechanics.  In  this  case,  for 
a  spherical  indenter  in  contact  with  a  sample,  the  radius  of 
the  contact  area,  a,  is 

l3FR\m 


where  F  is  the  applied  load  and  R  is  the  radius  of  curvature 
of  the  tip.  Equation  (6)  is  valid  for  a  nonsliding  contact.  For 
sliding  contacts,  Johnson45'46  predicted  that  the  contact  ra¬ 
dius  may  be  reduced  compared  to  static  contact.  However, 
the  experimental  evidence  using  AFM  by  Carpick  et  al.41 
and  surface  force  apparatus  (SFA)  by  Israelachvili48  suggests 
that  there  is  only  very  little  or  no  change  in  contact  area 
during  sliding  in  AFM  and  SFA  experiments,  respectively. 
Assuming  that  sliding  during  scanning  does  not  significantly 
change  the  contact  area,  Eqs.  (4)  and  (6)  can  be  combined  to 
obtain  the  reduced  modulus 


E*  = 


(7) 


Figures  6(a)  and  7(a)  show  the  reduced  modulus  images 
(storage  and  loss)  calculated  by  substituting  the  contact  force 
and  stiffness  data  into  Eq.  (7).  The  contrast  clearly  shows 
that  the  modulus  at  the  center  of  the  fiber  is  less  than  at  the 
periphery,  while  the  modulus  of  the  epoxy  is  at  least  an  order 
of  magnitude  less  than  the  fiber.  Edge  artifacts  can  be  clearly 
distinguished  from  true  mechanical  response,  but  are  mini¬ 
mal  here  due  to  a  good  sample  polish.  To  reduce  the  noise, 
the  data  were  averaged.  A  cross  section  of  the  data  running 
across  the  center  of  the  fiber  is  shown  in  Figs.  6(b)  and  7(b). 
The  storage  modulus  of  the  epoxy  is  —1.3  GPa,  while  the 
fiber  center  is  —35  GPa,  increasing  slowly  to  a  constant 
value  of  —53  GPa  at  the  periphery.  Loss  moduli  for  both 
regions  are  about  1/5  to  1/3  that  of  the  storage  moduli. 


C.  Comparison  with  standard  techniques 

For  comparison,  conventional  indentation  experiments 
were  done  at  different  locations  of  the  sample.  Figure  8 
shows  the  indentation  response  at  the  center  (A)  and  periph¬ 
ery  (B)  of  the  fiber,  and  in  the  epoxy  (C)  for  a  maximum  load 


of  200  /xN.  The  modulus  calculated  from  load-displacement 
graphs  using  standard  unloading  slope  analyses  is  3.5 ±0.5 
GPa  for  the  epoxy  region,  48 ±2  GPa  for  the  center  of  the 
fiber,  and  57  ±3  GPa  at  the  periphery.  This  confirms  a  —20% 
difference  in  modulus  across  the  fiber,  which  is  clearly  vis¬ 
ible  in  the  images. 

A  second  set  of  indentation  experiments  were  done  to 
compare  dynamic  properties  of  the  sample.  In  these  experi¬ 
ments,  storage  and  loss  moduli  were  calculated  from  a  con¬ 
tinuous  measurement  of  contact  stiffness  at  a  given  fre¬ 
quency.  Figure  9  shows  the  contact  stiffness  response,  at  a 
frequency  of  105  Hz,  as  a  function  of  load  for  the  carbon 
fiber  and  epoxy  matrix.  Because  during  these  experiments, 
the  indents  penetrate  the  sample  to  greater  depths  than  during 
the  contact  stiffness  imaging  experiment,  the  fit  for  the  Berk¬ 
ovich  tip  is  used  rather  than  the  Hertzian  approximation  to 
calculate  elastic  modulus. 

Figure  10  shows  a  comparison  of  modulus  measure¬ 
ments  by  indentation  and  imaging.  To  get  a  statistical  esti¬ 
mate  of  the  modulus  from  the  imaging  technique,  histograms 
of  the  image  data  were  analyzed,  and  mean  and  standard 
deviation  of  moduli  for  the  various  regions  were  obtained  by 
peak  fits  using  normal  distributions. 

The  reduced  moduli  calculated  from  the  scanned  images 
are  in  excellent  agreement  with  the  indentation  measure- 


Depth,  nm 


FIG.  8.  Conventional  load-displacement  indentation  curves  of  the  carbon- 
epoxy  composite  at  the  center  (A)  and  periphery  (B)  of  the  fiber  and  in  the 
epoxy  (C)  matrix  are  shown.  A  maximum  load  of  200  /rN  was  applied. 
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FIG.  9.  Plots  of  the  storage  and  loss  contact  stiffness  (in  N/m)  for  (a)  the 
epoxy  matrix  and  (b)  carbon  fiber  as  a  function  of  load  (in  /x N)  are  shown. 


ments  for  both  the  storage  and  loss  components.  Some  varia¬ 
tion  of  modulus  between  the  new  image  mode  and  conven¬ 
tional  indentation  mode  maybe  attributed  to  differences  in 
the  frequency  response  of  the  materials;  in  any  case,  the 
agreement  is  quite  good.  The  larger  scatter  in  the  moduli 
obtained  from  the  fiber  region  with  the  imaging  technique  is 
due  to  noise,  which  can  be  reduced  by  increasing  the  ampli¬ 
tude  of  the  force  modulation.  However,  imaging  conditions 
were  optimized  to  get  a  measurable  displacement  response 
from  both  fiber  and  epoxy  regions.  The  storage  moduli  com¬ 
pare  favorably  to  literature  values26  for  other  carbon  fibers. 
The  good  correlation  between  the  imaging  (low  load)  and 
indentation  (higher  load)  modes  suggests  that  the  sample 
polishing  technique  did  not  damage  the  surface,  as  was  ob¬ 
served  in  previous  low  load  indentation  experiments  on  ara- 
mid  fibers.49 


FIG.  10.  Comparison  of  modulus  measurements  made  by  indentation  and 
imaging  at  the  center,  (A),  and  periphery,  (B),  of  the  fiber  and  epoxy,  (C)  are 
shown. 


The  imaging  results  can  improve  the  way  we  conduct 
mechanical  property  measurements  in  several  ways.  Rather 
than  depending  on  topography  alone  in  a  scanned  image  to 
select  indentation  locations,  the  contrast  in  the  phase  and 
amplitude  response  can  be  used  to  identify  regions  of  inter¬ 
est.  This  could  be  especially  valuable  in  nanocomposites, 
polymers  or  tribological  surfaces,  for  example,  where  topog¬ 
raphy  may  not  correspond  to  mechanical  property  differ¬ 
ences,  or  in  studies  of  subsurface  delamination  of  thin  films. 
At  the  scan  rates  used  in  these  experiments,  the  images  take 
about  2-3  h  to  complete.  While  this  is  considerably  longer 
than  conventional  scanning  rates  in  AFM,  the  images  can  be 
analyzed  and  converted  to  quantitative  stiffness  and  modulus 
maps  in  a  matter  of  minutes.  Importantly,  this  method  elimi¬ 
nates  the  tedious  task  of  performing  and  analyzing  indenta¬ 
tion  arrays  where  sample  nonuniformities  (e.g.,  defects,  as¬ 
perities,  and  edge  effects)  are  difficult  to  locate  and  interpret. 
Furthermore,  because  the  tip-sample  interaction  is  elastic, 
the  spatial  resolution  is  also  improved  over  indentation  maps, 
where  indents  should  be  spaced  far  apart  to  eliminate  the 
effects  of  residual  stress  and  plastic  deformation  from  neigh¬ 
boring  indents.  Any  inelastic  tip-sample  interactions  (i.e., 
plastic  deformation)  will  also  limit  the  application  of  this 
technique;  hence,  experiments  should  be  performed  at  the 
lowest  practical  load.  The  results  confirm  that  quantitative 
mapping  of  dynamic  mechanical  properties,  with  the  tip  in 
continuous  contact  with  the  surface,  is  possible  if  proper  in¬ 
strument  calibration  and  analysis  are  performed. 

IV.  CONCLUSIONS 

Quantitative  imaging  of  contact  stiffness  and  elastic 
modulus  (storage  and  loss)  was  demonstrated  for  a  fiber- 
epoxy  sample  using  a  hybrid  nanoindentation  apparatus  that 
combines  depth-sensing  indentation  with  the  imaging  capa¬ 
bilities  of  AFM.  Good  correlation  of  storage  and  loss  moduli 
was  found  between  imaging  and  measurements  made  using 
conventional  nanoindentation  techniques.  Contrast  in  the 
phase  and  amplitude  images,  without  the  detailed  processing 
employed  in  this  article,  can  also  be  used  to  identify  regions 
of  interest  for  indentation  measurements  as  well  as  aid  in 
interpreting  subtle  changes  in  mechanical  properties. 
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APPENDIX 

The  transducer  and  its  accompanying  electronics  have 
been  designed  to  achieve  a  balance  between  maximizing 
signal-to-noise  ratio  and  preserving  bandwidth  sufficient  for 
quasi-static  operation  over  a  meaningful  range  of  loading 
and  displacement  rates.  In  the  case  of  dynamic  operation, 
however,  much  more  detailed  attention  must  be  given  to  fil¬ 
ters  (component  transfer  functions),  as  they  affect  the  mea¬ 
sured  response  versus  frequency  spectrum.  Of  course,  not  all 
transfer  functions  need  be  considered  in  terms  of  frequency 
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FIG.  11.  Block  diagram  of  relevant  filters  affecting  the  phase  and  amplitude 
of  the  ac  displacement  signal  is  shown. 


dependency,  rather  only  those  that  ultimately  impart  an  un¬ 
acceptable  level  of  error  to  the  determination  of  materials 
properties  over  the  dynamic  testing  bandwidth. 

Three  low-pass  filters  relevant  to  the  dynamic  testing 
bandwidth  used  in  this  work  exist  along  the  signal  path  from 
the  point  of  force  actuation  to  the  point  of  measuring  the 
displacements  (see  Fig.  11).  They  are  (1)  a  first-order  filter 
affecting  the  oscillatory  component  of  the  force  actuation 
voltage,  (2)  a  first-order  filter  after  the  synchronous  demodu¬ 
lator  that  is  integral  to  the  capacitive  sensing  technique,  and 
(3)  a  user-selectable  second-order  filter  operating  on  the  de¬ 
modulated  and  filtered  displacement  signal  just  prior  to  digi¬ 
tization  by  the  data  acquisition  system.  Transfer  functions  for 
these  filters  are  described  below  in  terms  of  RC  time  con¬ 
stants,  resistor  ratios,  and  frequency  a>  in  order  to  illustrate 
their  functional  form.  Specific  component  values  are  not 
given  as  they  have  evolved  over  the  course  of  time,  and 
would  not  be  universally  applicable. 
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with  the  phase  lag  introduced  to  the  signal  as 
tan  6=  coRC. 


(Al) 


(A2) 


Interestingly,  application  of  the  square-root  circuit  elimi¬ 
nates  coupling  between  ac  and  dc  components  in  determin¬ 
ing  the  force  amplitude  of  the  u>  term,  which  would  not  be 
the  case,  in  general,  for  arbitrary  a0  and  a , .  It  is  noted, 
however,  that  the  square-rooted  functional  form  deviates  sig¬ 
nificantly  from  a  simple  dc-shifted  sinusoid  in  the  limit  of 
Vac  approaching  Vdc,  therefore,  the  harmonic  content  must 
be  considerable  in  this  case.  Consequently,  knowing  the  re¬ 
sultant  coefficients  of  the  force  amplitude  spectrum,  in  the 
presence  of  the  actuation  filter,  takes  on  a  measure  of  com¬ 
plexity  in  that  limit. 

2.  Correcting  the  displacement  signal  (step  1) 

The  first  correction  to  the  measured  displacement  X  in¬ 
volves  the  transfer  function  of  the  active  first-order  low-pass 
filter  acting  on  the  output  of  the  synchronous  demodulator, 
which  is  described  by 


X, 


d—  output 


mechanical 


1+Rt/R, 


V(l+/?2//?1)2+(2^2C)2’ 


(A3) 


and  a  phase  lag  of 


1.  Correcting  the  applied  force 

The  force  modulation  electronics  used  in  this  work  em¬ 
ploy  a  square-root  circuit  to  generate  the  voltage  to  be  am¬ 
plified  and  filtered  by  the  actuation  electronics,  with  its  out¬ 
put  given  by  VSR=  'JVdc+  Vac  sin  cot.  This  circuit  is 
appealing  from  the  standpoint  that  for  a  parallel-plate  capaci¬ 
tor  the  electrostatic  force  is  proportional  to  the  square  of 
voltage,  which  would  lead  to  the  “clean”  result  of  the  ap¬ 
plied  dynamic  force  always  proportional  to 
Vdc+Vac  sin  cot.36  The  presence  of  the  actuation  filter,  how¬ 
ever,  precludes  straightforward  application  of  this  concept 
across  the  dynamic  testing  bandwidth. 

Typically,  the  transducer  is  biased  with  a  dc  voltage  Vdc , 
which  is  set  approximately  two  orders  of  magnitude  greater 
than  Vac;  hence,  evaluating  the  case  of  Vdc8>Vac  is  most 
appropriate  to  this  study.  In  this  limit,  the  square-rooted 
functional  form  (V^+V^  sin  wf  can  be  approximated  by 
a0  +  a  i  sin  wf,  where  a0  =  0.5[Vvdc+ Vac+  \lvAc—Vac]  and 
a1  =  0.5[Vvdc+  Vac—  \/Vdc—  Vac]  if  one  imposes  equal  oscil¬ 
lation  amplitude.  For  the  values  of  Vdc,Vac  used  in  these 
experiments,  this  approximation  deviates  maximally  from 
the  true  functional  form  by  only  0.001%;  the  harmonic  con¬ 
tent  at  2u>  is  nearly  three  orders  of  magnitude  smaller  than 
the  fundamental.  The  practical  implication  of  this  exercise  is 
that  the  oscillatory  component  of  the  dynamic  force  is  lim¬ 
ited  to  the  fundamental,  and  that  the  frequency  dependency 
of  its  amplitude  is  governed  by  the  transfer  function  for  the 
classical  passive  first-order  low-pass  filter  acting  on  a  sinu¬ 
soidal  input.  In  other  words, 


tan  0= 


2  (oRiC 

1  +R2!R\ ' 


(A4) 


In  Eq.  (A3),  AjrK.chanical  is  the  actual  displacement  amplitude 
and  Xd_output  is  the  filtered  displacement  amplitude  after  de¬ 
modulation. 


3.  Correcting  the  displacement  signal  (step  2) 

The  second  and  final  correction  to  the  measured  dis¬ 
placement  involves  the  transfer  function  for  the  user- 
selectable  active  second-order  low -pass  filter  just  prior  to  the 
point  of  digitization,  which  is  described  by 


y 

^  output 

Xr  • 

—  input 


l 

\/a2+b2  ' 


(A5) 


where  X/_input=Xd_output  from  Eq.  (A3),  A  =  1 
—  arR\C \R2C2  and  B  =  a >C2(/?]  +2 R2),  with  the  phase  lag 
given  by 


tan 


(A6) 


In  summary,  the  force  amplitude  is  governed  by  Eq. 
(Al),  the  mechanical  displacement  amplitude  can  be  found 
by  dividing  the  measured  displacement  amplitude  by  Eqs. 
(A3)  and  (A5),  and  the  mechanical  phase  shift  corresponds  to 
the  phase  shift  measured  relative  to  the  square-root  circuit 
output  minus  the  sum  of  electronic  phase  shift  (phase  lag) 
given  by  Eqs.  (A2),  (A4),  and  (A6). 
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D.  Alternative  (simplified)  phase  shift  correction 
method 

An  alternative  method  to  correct  for  phase  shift  due  to 
electronics  can  be  accomplished  by  zeroing  the  phase  at  the 
lock-in  amplifier  for  a  given  (fixed)  driving  frequency.  In  this 
case,  it  is  necessary  to  know  the  dynamics  of  the  transducer 
when  it  is  free  of  contact.27  The  corrected  phase  shift  <bc  in 
this  case  is 

(A7) 

where  C,  is  the  damping  coefficient  of  the  transducer,  K,  is 
the  spring  stiffness,  m  is  the  mass  of  the  transducer,  and  6  is 
the  measured  phase  shift  from  the  lock-in  after  zeroing  the 
phase  manually  on  the  lock-in.  In  the  present  experiment,  the 
phase  shift  due  to  electronics  is  66°  at  a  driving  frequency  of 
200  Hz. 
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